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Abstract 

If a Hamiltonian dynamical system with n degrees of freedom admits m constants of motion more 
than 2n — 1, then there exist some functional relations between the constants of motion. Among these 
relations the number of functionally independent ones are s — m ~ (2n — 1). It is shown that for such a 
system in which the constants of motion constitute a polynomial algebra closing in Poisson bracket, the 
Nambu brackets can be written in terms of these s constraint functionals. The exemplification is very 
rich and several of them are analyzed in the text. 

1 Introduction 

The concept of generalized Hamiltonian dynamics arouse in 1973 with an article 
by Y. Nambu. [1] In his proposal, Nambu employed an A^-ary bracket, generically 
called Nambu bracket (NB), to describe the time evolution of the dynamical system 
in A^-dimensional (ND) phase space. His bracket includes — 1 functionally 
independent constants of motion, the so-called generalized Hamiltonians. As an 
illustrative example, Nambu considered the Euler equations of free rigid body for a 
3D phase space and this was the only example given. Finding examples in higher 
odd-dimensions is still very tedious matter. 

In Nambu formalism, dynamical systems produce inevitably a nontrivial normal- 
ization factor C at least when N is an even integer grater than three. [2] In words, 
in order to get the correct Hamiltonian dynamics the NBs must be normalized 
properly. The explicit general form of C has been derived in detail for superinte- 
grable systems. [3] The aim of this paper is to obtain C for a Hamiltonian system 
with n degrees of freedom and m constants of motion Ci, • • • , C^, with m > 2n. 

First, we begin by reviewing the basic features of the Nambu formalism. Let 
M be an ND smooth manifold and let C°°(M) be the linear space of smooth 
real-valued functions on M. The real multilinear map 

{, } : C°°(M) X • • • X C^{M) C°°{M) (1) 
" V ' 

N times 

defines NB of A-th order satisfying the properties skew-symmetry, Leibniz rule and 
generalized Jacobi identity (fundamental identity). [4] When one considers this NB 
structure C°°(M) admits another algebra structure. 
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The time evolution of / G C°°(M) is determined by A^ — 1 Hamiltonian functions 
Ci, . . . , Cn-1 G C°°(M) and is described by the Nambu-Hamilton (NH) equations 
of motion 

| = X^^(/) = {/,Ci,...,CAr_i}, (2) 

where is called the NH vector field corresponding to Ci, Cat- i- 

Finally, we recall that when M is a symplectic manifold of dimension N = 2n, 
C°°{M) has also an infinite-dimensional Lie algebra structure defined with respect 
to the Poisson bracket (PB) 

n 

{/, 9}p = E (d,Jd,,g - dpjd,^g) , (3) 

where (q, p) = (gi, . . . , pi, . . . , pn) are the local Darboux canonical coordi- 
nates. 



2 Canonical Nambu Bracket and its Decomposition 

A concrete realization of NB was embodied with the following form 

1/1, • • • , jN} - -W7 V 

d{xi, ...,xn) 

by Y. Nambu in the case of M = and it is called the canonical NB. In (4) 
X = {xi, . . . ,xn) denotes the local coordinates of and the right hand side 
stands for the Jacobian of the mapping / = (/i, . . . , f^) : R^ R^. 

For the systems with n degrees of freedom ,i.e., 2n-dimensional phase space, the 
canonical Nambu bracket is determined explicitly by 

{/i, • • • , /2n} = -^7--^ ^T^^TT ^ ^h-iin^hh ■ ■ ■ di^nhn (5) 

and it can be decomposed into a skew-symmetric product of PBs which is a useful 
identity 

{/l, • • • , f2n} = {fu /[2}p{/3, • • • , f2n]}, (6) 

where the bracket [ ] indicates the cyclic sum with respect to its elements. [5] (Here 
and hereafter we will use the Einstein summation convention for all repeated in- 
dices.) 

For coordinate-free expression of the canonical NB we associate the (A/" — l)-form 

F = dCi A • • • A dC7v-i (7) 
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to given N — 1 Hamiltonian functions Cj. In (7), d and A denote the usual exterior 
derivative and exterior product of Cartan calculus. If we employ the Hodge map 
* : Ap(RN) AN-p(RN) between the space of j:^forms A^(R^), p < N and the 
space of {N — p)-forms A^~^(R^), an easy calculation shows that 

*(d/ A r) = {/. Ci. . . . , Cm-i} = ^^y^' • • • ' . (8) 

Finally, for the purpose of this study, we write the NH equation for any / as 

*{df A r) = C/, (9) 

where C is also a constant of motion. C will be referred to as normalization constant 
corresponding to the set Ci, . . . , Cn-i and will be specified from the requirement 
that (9) produces the correct Hamiltonian equations of motion in the case of = 
2n. The requirement of nontrivial C is inevitable at least when n > 2. The NB 
given in (9) will be used in the rest of the text. 



3 Nambu Brackets and Constraint Functionals 

While some authors have investigated the connection of Nambu dynamics to Dirac's 
constraint formalism, [6, 7, 8] our approach will depend on functions of these con- 
straints. We consider a system with n degrees of freedom, i.e., a 2n-dimensional 
phase space with coordinates (gi, . . . , n > 2. We suppose that the 

system acquires m> 2n constants of motion Ci, . . . , Cm constituting a polynomial 
algebra closing in PB. Though these restrictions seem a bit rigorous, it may be 
remarkable to emphasize that most of the examples in the literature are of this 
kind. We impose no restriction on involutive properties of the constants of motion. 
In this sense, the system need not be integrable or superintegrable. For such a sys- 
tem, when possible, the number of functionally independent constants of motion is 
2n — 1 and therefore there are independent m — {2n — 1) = s functional relations 

= FjiCi, ...,Cm)=0, (i = 1, . . . , s) (10) 

between C^s [9]. We will call the FjS constraint functionals to avoid any confusion. 
In some cases, too many constraint functionals may appear, but any independent 
s of them are enough to operate the formalism, hence the choice of the constraints 
for these systems is not unique. We divide the derivation into two main cases to 
obtain a full treatment. First, for the case s < 2n — 1, let us choose the independent 
set of Fi , . . . , Fg and construct the s-form 

3F 3F 

a = dFiA---AdFs = ^---^dCi,A---A dQ^. (11) 
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Now we take the arbitrary set Ci^_^^, . . . , Ci^^_^ and construct the ((2n — 1) — s)-form 

(3 = da^^, A ■ ■ ■ A da,^_,, (12) 

and then the (2n — l)-form 

a AP = dFi A - ■ ■ A dFg A dCi^^^ A - ■■ A dCi^^_^ 
dF dF 

= ■ ■ ~ dCi, A ■ ■ ■ A dCi^ A dCi^^, A ■ ■ ■ A dCi,^_,. (13) 

When we multiply (13) by df and apply the Hodge map, we get a 2n-th order NB 

^{df A a A /?) = {/, Fi, ... , F„ C,,,,, . . . , J. (14) 

Obviously, if the FjS are written in terms of the phase space coordinates then they 
get zero, therefore 

3F 3F 

*idf Aa A = • • • -^^idf AdQ, A ■ ■ ■ AdQ^^J = 0. (15) 

Each NB in 15, by the definition (9), implies a normalization constant associated 
withCi,,...,C,2„_i, i.e., 

'{df AdQ,A---A rfQ,„_J = fCi,..i,^_,. (16) 

Thus this leads to a linear homogeneous system 

dFi dF, 



da, oa^ 



Cn-i.n-. = (17) 



(TTl \ 
2n I si ^^^y basic solution techniques 

that (17) has infinitely many solutions for the unknowns Ci,:.i2n-i- -^^^ require 
a solution which is nontrivial and compatible with the decomposition (6). The 
suitable solution can be chosen as 

_ m_ dF, d(Fr,...,F,) 

where the sign rh is determined by the Levi-Civita tensor in the second term. 
Indeed, this choice argues with (17), 

dFi OF, dFi OF, 
^iv-isis+i ■ ■ ■ i2n-ii2n-im fjfi. ' ' ' f)n. pjn. f^n. 

(2n— 1 — s) indices 

^ (_l).(2n-i-.)^. OF, OF, 



= ^ a(Fi,...,F„Fi,...,F,) _^ ^^^^ 
d{Ci,, . . . , Ci^, Ci^^, . . . , Ci^) 
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An illustrative example may be convenient to be more explicit. Consider the system 
with two degrees of freedom including the constants of motion Ci, . . . , C5, which 
is also considered in Subsec.4.1.. Thus there exist two independent constraint 
functionals, Fi, F2. Therefore 



dFi dF2 
dC~dC~ 

diFuF2] 



a = dFiA dF2 = dd, A dQ^ 



dCi, A dCi,, {ii < 22); (20) 

and f3 = dCi^. For an arbitrary Ci^, say Ci, the condition *((i/A(iFi A(iF2A(iCi) = 
implies 

diF^ d{F,,F2) d{F,,F2) 

d{C2, C3) diC2, C4) diC2, C5) 

^ diF^,F2) d{F,,F2) d{F,,F2) 

Referring back to (18), one computes easily 

d{Fi,F2) d{Fi,F2) _ d{F\^d{F\^ d{F\^d{F\^ 

5(C2, C3) a(C4, C5) 5(C2, C4) 5(C3, C5) 5(C2, C5) 5(C3, C4) 

a(Fi,F2,Fi,F2) 



5(C2, C3, C4, C5) 



= 0. (22) 



At the second stage, i.e., for the case s > 2n — 1, we choose any 2n — 1 constraints 
from the set Fi, . . . , F^. With the same argument followed in the previous case, it 
is easy to construct the 2n-form 

\df A dF,^ A • • • A rfF,,„_J 

= |J---g^*(#ArfC,A...A^C,„J=0 (23) 
generating the system with I r> ^ 1 1 equations 







^ ( 2n - 


..) 


dF,, 









= 0. (24) 

Since the choice of the set Fj^, . . . , Fj^ _^ is arbitrary, (18) is always a solution to 
(24): 

dF^ dFj,^_, OF, dFs 
, , . . . , F- , Fi, . . . , Fo) 
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Consequently, if (16) is recombined with (18), one concludes that 
'^(d/ A dQ, A • • • A tiQ,„_J = {/, a„ . . . , J 

^il'"*2n-l*2n'"*m Q(J . OC ' 

*2n Im 

^ ^ d{Fi,...,Fs) j ^^^^ 

The formalism (26) possesses several NBs all are in accordance with correct equa- 
tions of motion. It may be remarkable to emphasize that we have no any restriction 
about the independence of the set (7^2^^, . . . , Ci„^. Thus, even if they are not inde- 
pendent, this does not destroy the validity of the formalism. On the other hand, 
the general result (26) also justifies the statement: If a NB includes any dependent 
subset of the constants of motion, then it vanishes. We shall make these remarks 
more precise in the discussion of the examples. 

Now, as a corollary, we conclude the statement: If one of the constants of motion, 
say Cfc, is taken as the Hamiltonian, then the decomposition (6) can, by discarding 
the Hamiltonian, be written in terms of the constraint functionals. The proof is 
straightforward: By the decomposition (6), 

{/, Ck, Ci^, . . . , Ci^^_^} = {/, Ck}p{Ci^, . . . , Ci2„_2} = f{Ci^, . . . , Ci^^_^}, (27) 

on the other hand, by (26), 

dFi dFg ■ 

{/? C'fc, Cj^, . . . , Cj2„_2} = ^kiv-im-i 'f)r^, ' ' ' ^ ' ^^^-^ 

*2n— 1 I'm—I 

thus 

5(Fi, . . . , Fg) 

{Ci^, . . . ,Ci^^_^} = — r. (29) 

In particular, for the case n = 2, (29) holds for the PBs. 

Finally, we talk about determination of the FjS. Note that, as a general result 
of (6), for the NBs including the Hamiltonian, the normalization constant C is 
obtained easily in terms of the PBs. Therefore this observation supplies a helpful 
guide in determining the constraint functionals. Under this circumstance, s — 1 
functional relations can be taken as the constraint functionals without any rear- 
rangement. Thus, the construction of the last one is reduced to the problem of 
finding a function whose derivatives are known. Since the PBs of the constants of 
motion close to constitute a polynomial algebra, this gives us an easy integration 
process. 
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Having shown how to construct such a formaUsm we will give examples to confirm 
its correctness. To be more clear, the first example (harmonic oscillator) has been 
analyzed in detail. The examples have been chosen in a variety so that they include 
various alternatives when considering the dimension of the phase space and the 
number of the constraint functionals. In all examples, although the construction 
of constants of motion is not unique, we kept the forms and numbers of them just 
as appeared in the literature cited in the text. 



4 Systems with Two and Three Degrees of Freedom 
4.1 Hcirmonic oscillator 

Although we take the system as the one with two degrees of freedom, its 2n, (n > 2) 
dimensional extension is obtainable as a special case from the Winternitz system 
given in Subsec.5.1. The system is described by the Hamiltonian 

Ci = H = pV2 + A;qV2, (30) 

where = gi^ + g| and = pi^ + pi [2]. Suppose that in addition to the 
Hamiltonian we are given the following set of equations as the constants of motion, 

C2 = pl/2 + kql/2, Cs=pl/2 + kql/2, 

Ca = qiP2-q2Pu = PiP2 + kqiQ2- (31) 

Their nonvanishing PBs are given by 

{C2,C4}p = — {C3, C4}p = — C5, 

{C2,C^}p = —{C's,C5}p = kC4, 

{CC^jp = -2(C2-C3). (32) 

If the functional relation Ci = C2 + C3 is taken as the first constraint such that 
Fi = Ci — C2 — C3, then the other which is compatible with the all PBs via (29) 
can be constructed as the following 

F2 = 2C2C3 - ^kCl - hol (33) 



so that 



r\ - ^(^• ^) _dF2_ 
d{Fi F2) dF2 
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Now we list two of the NBs as the sample, 



d{Fi F2) ■ 

{/,Ci,C2,C3} = |U3y/ = o (35) 

which are consistent with 

{/, Ci, C2, C4} = {/, Ci}p{C2, C4}p = —C^f, 

{/,Ci,C2,C3} = {/, Ci}p{C2, Csjp = 0. (36) 

On the other hand, for the NBs not including the Hamiltonian, C is not so evident 
via (6). For example 

{/, C3, C4, C5} = {/, Cs}p{C,, C,}p + {/, C5}p{C3, C,}p + {/, C4}p{C5, Csjp 

= -2(C2 - Cs){f, Cs}p + C5{/, Csjp + kC,{f, C,}p, (37) 
but it is straightforward to obtain it by the virtue of 

{/, C3, C4, C,} = = 2C3/. (38) 

If we turn off the Hamiltonian, the set C2, C3, C4, C5 is still closed in PB and we 
need only one constraint functional, F = —F2, thus 

{/, Cs, C4, C5} = = 2C3/. (39) 

Now consider the bracket 

{/, C^, C4, C5} = = -2(C2 - C3)/, (40) 

and take artificially Cq = C2 — C3 so that C2, C3 and Ce are not independent. Again 
{/. Ci. C4, C5} = f J'^. / = -2iC, - Of, (41) 

0*1^^6, ^2, ^3) 

where F3 = Ce — C2 + C3. 

Finally, for the case s > 2n — 1, let C7 = C4C5. Of much NBs, we have chosen 

where F4 = C7 — C4C5. 
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4.2 Smorodinsky-Winternitz system 

Smorodinsky-Winternitz system consists of a set of four Hamiltonians which have 
potential form, i.e., H = + V [10]. All potentials are separable into at least 
two coordinate systems and they also admit superintegrable structure. We will 
consider symbolically only the system 

Ci = H = pV2 + u;2(4gi2 + q^^) + aiqi + a2/g2^ (43) 

where all Greek letters are some real constants, q = {qi-\- ^i)^'^^ = pi' +P2- 

All constants of motions are at most quadratic in momenta, 

C2 = + aigi + 4u;V^ 

C3 = 2p2Lz - ^(^^qiql + ^a2qi/ql - aiql, 

C4 = -2{ai + 8uj\i)q2P2-pii2pl-4:Uj\l + 4a2/ql), (44) 

where L3 = qiP2 — is the third component of angular momentum. Their 
nonvanishing PBs close in a Poisson algebra 

{C2, C'sjp = C'4, 

{C2,C4}p = 4aiC2 - Scj^Ca - 4aiCi, 

{C3,C4}p = -48C| + 64C1C2 - 4aiC3 + 64a;2a2 - 16Cl^ (45) 

admitting a Casimir 

F = C|/2 - 4aiC2C3 + 4a;2c| + 4aiCiC3 - 16C| 

+ 32CiC| + 6Wa2C2 - lQCfC2 + 4a?a2 (46) 

as the constraint functional [11]. Note that the constraint (46) can be obtained as 
the integration of the PBs in (45). This observation justifies the expressions 

OF OF dF 

{C2,C3}P=|^, {C2,C4}P = -^, {^3,C4}P = ^. (47) 

And hence two of possible NBs are 

dF ■ 

{/,Ci,C2,C4} = --— /=-(4aiC2-8a;2C3-4aiCi)/, 

dF ■ 

{/,C2,C3,C4} = -^/=-(4aiC3 + 32C|-32CiC2)/. (48) 

4.3 Kepler- Coulomb system 



Let us concentrate now on the 6D Kepler-Coulomb Hamiltonian 

H = pV2 - a/q, 



(49) 
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where a is a real constant, = pi^ + + q = {qi^ + + Qs^Y^'^ [2]- 

Because of the rotational symmetry of the system, the angular momentum L is 
integral invariant and hence its components can be taken as constants of motion, 
i.e., 

Li = q2P3 - q3P2, L2 = qsPi - qiPs, = qm - q2Pi- (50) 

(From now on, we will write the constants of motion just as appeared in the system 
without corresponding any Ci to them to avoid any confusion in mind and to keep 
their symmetries in writting). Moreover, there also exists an extra invariant arising 
from that the particle has a closed orbit. This invariant is called the Runge-Lenz 
vector A given by 

A = p X L — aq/g. (51) 
In addition to the previous constants of motion, 

Ai = P2L2, - P3L2 - aqi/ g, 
M = PsLi - piLs - aq2/ g, 

A3 = P1L2 - P2L1 - aq^/q. (52) 
Consequently we have seven constants of motion satisfying the commutations 

{La, Lh}p = ^abcLc, {^a, A}p = —'^HCabcLc, {La, Ab}p = eabc^c, (53) 

and the following functional relations 

A • L = 0, = 1 + 2HL'^ (54) 

which are candidates for the constraint functionals. Thus if we choose these func- 
tions as follows 

Fi = AiLi + A2L2 + AsLs, 

F2 = l^HiLl^Ll^LD-^iAl^Al^Al), (55) 
(29) implies several brackets such as 

{Li, L2, L3, A,} = = A2LS - A,L2. (56) 

This is compatible with the decomposition 

{Li,L2,Ls,Ai} = {Li,L2}p{L3,Ai}p^{Li,Ai}p{L2,Ls}p 

+ {Li,Ls}p{Ai,L2}p. (57) 

So, one of the NBs is 

{f,H,L,,L2,Ls,A,} = {f,H}p{L^,L2,Ls,A^} 

d{Fi,F2) 



diA2, A,) 



f = (A2L3 - A3L2)/. (58) 
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Other two of them may be chosen as the following 

{/,Li,L2,L3,Ai,A2} = ^i^il^/ = L3( V + + Ls V, 
{f,H,L^,L2,A2,A,} = Y''^^\ f = -{A,A, + 2HL^L^)f. (59) 



5 Systems with n Degrees of Freedom 
5.1 Winternitz system 

Winternitz system is the arbitrary dimensional generalization of one of the Smorodin- 
sky - Winternitz Hamiltonians mentioned above. Their constants of motion are 
constructed by using the Lax matrix representation [12]. The particle's Hamilto- 
nian in n degrees of freedom is given by 



1 n ^2 



H = ^Up! + k'x^ + '-^), (60) 



where k and the ki are real constants. (Throughout this section all subscripts range 
from 1 to n). The elements of degeneracy group SU{n) of the Winternitz system 
are taken as the constants of motion all commuting with the Hamiltonian. First 
group of these functions has the form 

Tii = ^{Hi-kki), (61) 
Ik 

where the conserved quantity Hi is taken as the energy in the i-th direction. Second 
group is given by the ansatz 

T,j = f{H,)f{H,)A,A;, i^j, (62) 

where 

2k \ 



and 



1 / 



A = 1^ I + ^ - ^ x\ + IxkxiVi 1 . (64) 



'I 



The PBs of the functions T^- argue 

{T^j, Tys\p = i6jrTis — \5isTj.j. (65) 
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With this argument, there exist totally + 1 constants of motion and there must 
be functional relations between the invariants. One of the functional relations, so 
by (61), is the simplest one 



n 



H = j:{2kTii + kki). (66) 

i=l 

After a series of calculations, the others can be expressed as 

TijTjk = TjjTik. (67) 

Note, by referring to the discussion in Sec. 3, that the number of independent con- 
straint functionals is — 2n + 2. 

After having been defined the Winternitz system and its invariants, we now 
perform the case of three degrees of freedom as an example. Despite the fact that 
it is possible to study with ten constants of motion and five constraint functionals 
(or ten functional relations), for the sake of simplicity, we prefer the set of constants 
of motion H, Tn, T22, Tg^, T12, T13 which is also closed in the PB. The suitable choice 
for the only one constraint functional is 



F — Ti2Ti^ 



+ Tn + + Tss + ]^{ki + k2 + ks) 



(68) 



Consequently among the all possible four nonvanishing NBs two of them are listed 
in the following 

dF ■ 

{f, Til, T22, T12, Tia} = -^=rf = TuTi^f, 

{/,rn,r22,r33,ri2,ri3} = -^f = jj^^uTisf. (69) 



5.2 Free peirticle on n-sphere 



Our last example is a free particle moving on the surface of an n-sphere with the 
radius q = (qi^ + • • • + qr^Y^'^ [5]. The PB Lie algebra of the integral invariants 
(charges of so(n + 1)) is generated by the n(n — l)/2 rotation elements L^/? = 
QaPp — qpPa, Ck;, = 1, . . . n, and the n momenta Pq, = (1 — q^Y^'^Pa- The PBs of 
the invariants are 

{Lq,/j, L^^]p = 6p^Laj + SajLp^ — djSjLa^ — da^LjSj, (70) 
{Pa, P/sjp = Lap, {La(3, P-f}p = ^ajPp " ^jSjPa- (71) 
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The Hamiltonian of the particle is given by 

H = ^{PaPa^L0^Lp^), ^<7. (72) 

Unhke the previous example, the algebra is closed by all of the n{n + l)/2 + 1 
constants of motion, and then we need n(n — 3)/2 + 2 independent functionals to 
proceed the formalism. For this aim we will use the following concluded functional 
relations as the source of the constraint functionals; 

-^[0102-^0:3-^04] 0? 

-^'[0102-^03-^04 ■ ■ ■ Pan] — 0- C^^) 

Additionally a second group appears as 

-^'oi [02 -^0304] = 0- C''^) 

For each of the relations in (73) and (74), the a^s can be chosen freely provided they 
are all different from one another. It is clear to see that the number of functional 
relations listed above exceeds the needed one too much, but any set of independent, 
suitable rearranged n{n — 3)/2 + 2 functionals does work. As was the case for the 
previous example, we restricted ourselves to a particular, say n = 4, degrees of 
freedom. In that case, we have 11 constants of motion and therefore we need four 
independent constraints. First one is the simplest one, i.e., (72), 

Fi = H - ^{Pi' + P2' + P,' + P4') 

— -(-^12^ + -^13^ + -^14^ + -^23^ + -^24^ + -^34^)- (75) 

For the other three, a suitable choice can be written explicitly as the following 

F2 = L12L34 + L14L23 — L13L24, 

^3 = ^12^3 — L1SP2 + ^23^5 

F, = P,-^P, + ^P,. (76) 

Now, as before, we choose a sample from the NBs, 

d{Fi, F2, F3, F4) ■ 
{/j Pi, P2, P3, Pa, L12, L13} = — — — — ' ' — r/ 

OyLilA, i^23? -^24, -i^34j 

= [^12^2^4 + Li^P3Pa - Lu{Pi^ + P2' + P3')] /, (77) 
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-P4' = t) ^ ^12 ~ PsLu + PiLsA H — ^^L23 (78) 



gives the same equations of motion. For the readers who may wonder of other 
brackets, here we hst two of them, 
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